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I 

We investigate the dynamics of a spherically symmetric vaccum on a Randall 
and Sundrum 3-brane world. Under certain natural conditions, the effective 
Einstein equations on the brane form a closed system for spherically symmetric 
dark radiation. We determine exact dynamical and inhomogeneous solutions, 
• which are shown to depend on the brane cosmological constant, on the dark 

, radiation tidal charge and on its initial energy configuration. We identify the 

conditions defining these solutions as singular or as globally regular. Finally, 
■ we discuss the confinement of gravity to the vicinity of the brane and show that 

, a phase transition to a regime where gravity is not bound to the brane may 

I occur at short distances during the collapse of positive dark energy density on 

' a realistic de Sitter brane. 
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1 Introduction 

It is possible to localize matter fields in a 3-brane world embedded in a higher di- 
mensional space if gravity is allowed to propagate away from the brane in the extra 
dimensions (see for example ^ and the early related work in P]). The extra dimen- 
sions were usually compactified to a finite volume but it was realized that they could 
also be infinite and non-compact as in the model of Randall and Sundrum (RS) who 
proposed gravity to be confined to the vicinity of the brane by the warp of a single 
extra dimension 0. 

In the RS brane world scenario the observable universe is a 3-brane boundary of 
a non-compact Z2 symmetric 5-dimensional Anti-de Sitter (AdS) space. The matter 
fields are restricted to the brane but gravity exists in the whole AdS bulk. The 
classical dynamics is defined by the 5-dimensional Einstein field equations with a 
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negative bulk cosmological constant and a Dirac delta source representing the brane. 
The original RS solution involves a non-factorizable bulk metric for which the light- 
cone into the fifth dimension is closed by an exponential warp factor B 



7.2 



ds 



where 7]^^ is the Minkowski metric in 4 dimensions, 7 = y—An'^/Q with A the negative 
bulk cosmological constant and k"^ = Svr/Mp with Mp denoting the fundamental 
5-dimensional Planck mass. The brane cosmological constant was assumed to be 
fine-tuned to zero and so A = — k^A^/G or 7 = where A is the brane tension. 

Weak field perturbations around this static vaccum solution 0-0 were also analysed 
to reveal a new way to solve the hierarchy problem and show that 4-dimensional 
Einstein gravity is effectively recovered on the brane at low energy scales if the AdS 
radius I/27 is small enough. In particular Newton's potential Vn generated by a 
mass M receives warp corrections which to leading order give 



V, 



GnM 



N 



1 + 



/-y2 ^2 



(2) 



where Gj^ is Newton's gravitational constant. Extensions of this scenario were rapidly 
developed to allow for example non-fine-tuned branes 0, H and thick branes . The 
Gauss-Codazzi formulation used in 



10 



^ was discussed earlier in 
Friedmann- Robertson- Walker cosmologies fill- 13, static black holes |ll5l-ll8| and 
stars [0, the Vaidya solution |^ as well as Oppenheimer-Snyder gravitational col- 
lapse solutions ||2l|, ^ were also discovered within the RS brane world scenario. 
However, it should be noted that the static black hole solutions discussed in are 
for a 4-dimensional bulk and the existence of black hole solutions in the bulk that 
reduce to the a static black hole localized on the brane remain unknown. Moreover, 
cosmological perturbations have started to be implemented p3[-|2^. So far a consis- 
tent adjustment to the available data coming out of astrophysics, cosmology and also 
high energy particle colision experiments has been shown to hold |2^. In addition, it 
has been shown that the AdS/CFT duality pSf is compatible with the RS scenario 

m. 



Studies of the gravitational collapse of matter have been limited to static or ho- 
mogeneous dynamical models. Notable examples are the 5-dimensional black string 



solution [|T5[, the tidal Reissner- Nordstrom black hole [|T3, the Vaidya solution on the 
brane [^] and the Oppenheimer-Snyder model j^, As a result it has become 
clear that the exterior vaccum of a collapsing distribution of matter on the brane 
cannot be static PTl This was shown to be a consequence of the presence of 
gravitational modes in the bulk, which act as a kind of intrinsically dynamical dark 
matter in the whole brane world. These gravitational degrees of freedom are generated 
by the existence of Weyl curvature in the bulk and are a source of inhomogeneities 
in matter clouds located on the brane. 
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The purpose of this paper is to analyse the spherically symmetric, inhomogeneous 
vacuum on the brane. This is a necessary step towards the understanding of a real- 
istic collapse problem in the RS scenario. We take an effective 4-dimensional point 
of view following the Gauss-Codazzi covariant geometric approach |10|, |ll|] (see 



also Ref. for a review and notation). In section 2 we introduce the Einstein 
vaccum field equations on the brane. Although these equations do not, in general, 
form a closed system, additional conditions may be imposed to form a closed system 
described by two free parameters, viz., the brane cosmological constant A and the 
dark radiation tidal charge Q. In this section we also present a discussion on the 
localization of gravity near the brane. For small but positive A we show that while 
for Q < the tidal acceleration in the off-brane direction is always negative, implying 
acceleration towards the brane, for Q > the acceleration changes sign at a critical 
short distance. In section 3 we present A and Q dependent inhomogeneous dynami- 
cal solutions. We show that they take the LeMaitre-Tolman-Bondi form and depend 
on the initial configuration of the dark radiation, which is parametrized by a single 
function, / = /(r) > — 1. This is interpreted as the energy function. The marginally 
bound solution (/(r) = 0) is the static, zero mass tidal Reissner-Nordstrom black 
hole solution of |T^. In section 4 we analyse the dynamics of the non-marginally 
bound, inhomogeneous solutions and characterize precisely how they depend on A, Q 
and f{r). We also identify the conditions under which the dark radiation dynamics 
leads to the formation of a singularity as the final outcome of collapse or to a globally 
regular evolution with a bounce developing inside the vaccum. We conclude in section 
5. 

2 Vaccum Field Equations on the Brane 

In the RS brane world scenario the Einstein field equations in the bulk are |27 



r.2 



Gab — 



-Aqab + S{y) {-XgAB + Tab) , (3) 

The 4-dimensional 3-brane is located aty = and is a fixed point of the Z2 symmetry. 
Its induced metric is qab = Qab — nA^B where ha is the spacelike unit normal to 
the brane. Matter confined to the brane is characterized by the energy-momentum 
tensor Tab and satisfies Tab^^ = 0. The 5-dimensional metric is of the RS form 

ds"^ = QABdx^dx^ = dy^ + g^^dx^dx^ , (4) 

where g^i, should display a dependence on the fifth dimension which localizes gravity 
near the brane and closes the light-cone. 

According to the effective geometric approach 0, [10|, ^ the induced Einstein 



field equations on the brane are obtained from Eq. (J^), the Gauss-Codazzi equations 
and the Israel conditions with Z2 symmetry. In the vaccum Tab is set to zero, as a 
consequence of which the induced equations take the form 
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G 



where 
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(5) 



(6) 



is the brane cosmological constant and £^i, is the hmit on the brane of the projected 
5-dimensional Weyl tensor, 



hm 5^5^ 8ab = hm 6^6^CACBDn''n 



It is a symmetric and traceless tensor due to the Weyl symmetries and is constrained 
by the conservation equations 



V 



0, 



(8) 



obtained from Eq. (^) as a result of the Bianchi identities. 

The system of vaccum field equations on the brane defined by Eqs. (|[) and (||) 
may naturally be interpreted as defining the interaction of 4-dimensional Einstein 



gravity with matter represented by the traceless energy-momentum tensor 



- ^iivi (9) 
where = Svr/Mp with Mp = Xj^/G^ the effective Planck mass on the brane and 



rA 
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Using general algebraic symmetry properties it is possible to write E^y as |25 



K 



(10) 



11^ 



where such that u^u^ 



— 1 is the 4- velocity field and h^y = g^y + u^Uy is the 
tensor which projects orthogonaly to u^. The forms U, V^y and represent the 
effects on the brane of the free gravitational field in the bulk. Thus, U is the effective 
energy density. 



U 



Vny is the anisotropic stress. 
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(12) 



(13) 
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and Qfj, is the effective energy ffux 

2m= (14) 

Even though S^j^^, is trace free and has to satisfy the conservation equations, in general 
it cannot be fully determined on the brane. This is to be expected because observers 
confined to the brane cannot predict all bulk effects without a knowledge of the 
solution of the full 5-dimensional Einstein equations. As a consequence, the effective 
4-dimensional geometric theory is not closed. However, under additional simplifying 
assumptions about the bulk degrees of freedom it is possible to close the system if S^^, 
is constrained in such a way that it may be completely determined by its symmetries 



and by the conservation equations. This is what happens when V^^, = |j25| or if 
Vfj^u 7^ when a spherically symmetric brane is static |T^. 

Let us now show that it is possible to take a non-static spherically symmetric 
brane with Vfj,u 7^ and still close the system of dynamical equations. Consider the 
general, spherically symmetric metric in comoving coordinates {t,r,6,(f)), 

ds^ = g^^dx'^dx'' = -e^df + A^dr'^ + R^dU^, (15) 

where dVt^ = dd"^ + sm^9d(j? , a = a{t, r), A = A{t, r), R = R{t, r) and R is interpreted 
as the physical spacetime radius. If there is no net energy flux, then = 0. Again, 
if the stress is isotropic, V^v will have the general form 

V^u = V [r^r, - ^ V) . (16) 

where V = V{t, r) and is the unit radial vector, given in the above metric by 
= (0, A, 0, 0). The projected Weyl tensor then takes the diagonal form 



K 



^1= - ) diag (p, -pr, -pt, -Vt) , (17) 



where the energy density and pressures are respectively p = U , pr = (1/3) {U + 2V) 
and pt = (1/3) {U — V). Substituting in the conservation Eq. (H) we obtain the 
following expanded system [^| 



2^(p + Pr) = -2p-4|(p + pr), 

R' 

^'(P + Pr) = -p/ + 4 — (PT-Pr), (18) 

where the dot and the prime denote, respectively, derivatives with respect to t and r. 
A synchronous solution is obtained by taking the equation of state p + = 0, giving 
V = —2U. Then the conservation Eqs. ([18|) separate and we obtain 
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U + &A = Q = U' + A^U. (19) 
R R 



Consequently the effective energy density has the dark radiation form [113 



where the dark radiation tidal charge Q is constant. Thus for the gravitational dy- 
namics of spherically symmetric but still inhomogeneous dark radiation £^i, is fully 
determined to be 

£,iv = -^{u^,Uy -2r^,ru + h^,y) . (21) 

Note that this projected Weyl tensor £^j^ is a simple dynamical extension of the 
corresponding form associated with the static tidal Reissner-Nordstrom solution on 



the brane . In spite of the existing pressures in the dark radiation setting we may 
safely take the synchronous comoving frame, for which a = 0. Then the intricate 
general field equations simplify and substituting Eq. (^) in Eq. (^ we obtain 



G^,y = -Ag^,u + iu^,u,y - '^r^.r^ + h^^) , (22) 

a consistent and exactly solvable closed system for the two unknown functions, A{t, r) 
and R(t, r). 

Clearly, the dark radiation dynamics depends on A and Q 2^- An important 



point to note is that these parameters have a direct influence on the localization of 
gravity in the vicinity of the brane (see also Ref. ||14[)- To understand how consider 
the tidal acceleration away from the brane as measured by brane observers ES]. For 



the dark radiation vaccum we are considering such acceleration is |jT7| 

- lim ^ABCDri^n^n^M^ = + (23) 
y~*o± b R^ 



where ua is the extension off the brane of the 4- velocity field which satisfies u^riA = 
and u^ua = — 1- For the gravitational field to be localized near the brane the tidal 
acceleration must be negative. The condition for this to happen is 

Ai?4 < _^ (24) 

Consequently, the localization of gravity for all R is only possible if A < and 
Q < or A = and Q < 0. In terms of brane parameters this implies that A < Ac 
with Ac = and Q < or A = Ac and Q < 0. For A < Ac and Q > 

the gravitational field will only remain localized near the brane if R > Rc where 
-Rc = 3(5/(Ac — A). On the other hand for A > Ac and Q < confinement is 
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restricted to the epochs R < Rc- If A > Ac and Q > then gravity is always free to 
propagate far away from the brane. 

Since the dark radiation term comes from the electric or Coulomb part of the 5- 
dimensional Weyl tensor the geometry of the bulk may assign to Q any real value. 
Note however that in our settings the weak, strong and dominant energy conditions 



3T| all independently imply Q > because p = U, pr = —U and pt = U- This is 



consistent with other studies of dark radiation effects [^] but not with perturbation 
theory and the negative induced energy condition [|18] which need Q < 



Current observations do not yet constrain the sign of Q ^ ^ So in what 
follows we will keep Q as real parameter. The same does not hold for the cosmological 
constant A. Indeed, according to the present data A ~ 10~*^GeV^ ||33|. On the other 



hand Mp > lO^GeV and Mp ~ lO^^GeV imply by Eq. (|T0|) A > lO^GeV^ [g6]. Since 
using once again Eq. (p!0D we may write A^ = then A^ has a lower limit, 

Ac > 10~^^GeV^. Hence according to observations A must be positive and below the 
critical value Ac, < A < Ac. This implies a bulk with negative cosmological constant, 
A < 0. Note as well that the same conclusion holds if Mp is in the TeV range as Ac 
increases when Mp decreases. Finally note that when < A < Ac and Q > then 
gravity is not confined near the brane for R < R^. This is a non-perturbative result 
which is consistent with the energy conditions |^ for dark radiation on the brane. 
So far experiments have not been able to detect any confinement phase transition to 
a high energy regime where gravity is not bound to the brane. If such a transition 
does occur at l/R^ ~ ITeV then for Ac ~ lO-^^GeV^ we get Q ~ IQ-^^GeV-^, a 
value well below current astrophysical or weak gravity experimental constraints on Q 
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For completeness in what follows we will also discuss the de Sitter (dS) models 
with A > Ac which correspond to a dS bulk, the brane AdS setting where A < and 
the brane with A = 0. 



3 Dynamical Inhomogeneous Solutions 

To completly define the brane world dark radiation vaccum dynamics we need to 
solve Eq. (p^ . The solutions determine the evolution of the brane world as a whole 
and so are cosmological in nature. In the proper frame cr = and so the remaining 
components of the spherically symmetric metric given in Eq. (|15|) must satisfy the 
off-diagonal equation 

Gtr = ^ {AR' - R'A) = 0, (25) 

which implies 

A = ^, (26) 
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where H = H{r) is an arbitrary positive function of r. Then introducing Eq. (|26|) in 
the trace equation 



-G\ + G: + 2Gl = {aR + 2RA) = -2A + ^ (27) 

and integrating in r leads to 

fl^fl+l- (28) 

Note that at this point we have chosen to set the arbitrary integration function (of 
time) to zero. We make this choice here because such a function would correspond 
to an initial dust mass. A further integration gives 

R" = ^R"-% + /, (29) 

where / = /(r) is an arbitrary function of r. This function is naturally interpreted 
as the energy inside a shell labelled by r in the dark radiation vaccum. Integrating 
Eq. (HD we obtain 

r RdR , , 

J ^iR^ + fR^-Q' 

where i{j = ijj{r) is another arbitrary function of r and the signs + or — refer respec- 
tively to expansion or collapse. 

As it stands the dark radiation is characterized by three arbitrary functions of r, 
namely, H, f and ip. However, we have the freedom to rescale the radial coordinate 
r. Then we may impose on the initial hypersurface t = the condition 

R{0,r)=r (31) 

and so prescribe ip to be given by 

^ = / , (32) 

where the r.h.s is evaluated at t = 0. Note that condition (^) also defines the initial 
dark radiation effective density profile. Using the initial distribution of velocities we 
can also determine the energy function /. is then determined by the remaining 
Einstein equation, 

G: = -a - (33) 



as H = y/1 + / and consequently the metric takes the LeMaitre-Tolman-Bondi form 
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-dr + 



R 



i2 



1 + / 



dr^ + R^dVL^ 



(34) 



where / > —1. 

Note that for the marginally bound models corresponding to / = this metric 
describes a static solution. Indeed performing a transformation from the Tolman- 
Bondi coordinates (t, r) to the curvature coordinates (T, R) such that 



we find 



T = t+ dR 



R^-Q + f i?^ 
^R^-R^-Q' 



(35) 



ds' 



1 + ^ 

i?2 



3 



+ R'^dn^, (36) 

which is the inhomogeneous static exterior of a collapsing sphere of homogeneous dark 



fi2 
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radiation ElU (see also Ref. 22 



When A = this corresponds to the zero mass 
limit of the tidal Reissner- Nordstrom black hole solution on the brane |]T^ . The single 
black hole horizon only exists for Q < and is located at Rh = \/ —Q- Consequently 
for Q > the singularity at -R^ = is naked. For non-zero A the solution has horizons 

P 



at 



R 



± 




(37) 



If A > and Q < then we have an inner horizon R^ and an outer horizon R\. The 
two horizons merge for Q = — 3/(4A) and for Q < — 3/(4A) the singularity a.t Rg = 
becomes naked. For A > and Q > there is a single horizon at R^. If A < and 
Q > the singularity at i^^ = is again naked. For A < and Q < this no longer 
happens as a single horizon forms at R^. 

In the absence of dark radiation we simply find the homogeneous dS or AdS spaces 
|3^ . In particular for A > the solution takes the Friedmann- Robert son- Walker 
steady state form where the spatial curvature k is zero, 



R = r exp 




(38) 



The epoch i? = is regular but 



Then it represents (half) homogeneous dS space , 
fixed as it is impossible to leave it classically. 

For A = the marginally bound models allow us to calculate the leading order 



contribution to Newton's potential. Following Ref. IjT^ and keeping the coordinates 
T and R we may introduce a mass M for the gravitational source to find the tidal 
Reissner-Nordstrom metric 
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^^ + -|) dR' + RHn\ (39) 



Then Newton's potential is given by 



V, 



GnM 



Q 



Contrary to Eq. (H) this is a short distance perturbative expansion [0. So just as 
in the static hmit our djTiamical solutions with non-zero / should be more adequate 
to the high energy regime as long as quantum gravity effects can be neglected. The 
potential shows that perturbatively it should be Q < if anti-gravity effects are to 
be avoided at short distances (see also Ref. [|l8l). Consequently if a positive Q is to 
be allowed as is by the standard energy conditions |^ than it must be seen as an 
intrinsically non-perturbative effect. 

Let us now consider / 7^ 0. Depending on / and on A and Q we distinguish several 
possible dynamical inhomogeneous solutions. 

3.1 dS dynamics: A > 

In this setting Eq. (pO| ) can be cast in the form 



±t + ^ = 7;\T 



dY 



2VAV y^Y-'-P' 
where Y = + 3//(2A) and /? = /?(r) is given by 



(41) 



P 



A 



3p 
4A 



Q 



(42) 



Direct evaluation of this integral determines R and shows that the solutions are 
organized by p. For /3 > we find 



^ +2A 



/3 cosh 



±2W + cosh"^ 



r2 _L 3/ 



VP 



(43) 



where for Q > the energy function may span all its range, / > — 1, but for Q < 



it must further satisfy |/| > 2J—QA/3. If /? < we obtain 



2A 



-f3 sinh 



±2W^t + sinh"^ 



V2 -U 3/ ^ 
' 2A 



-P 



(44) 



where Q < and |/| < 2J-QA/3. For /3 = we get 
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^2 3/ 



2^3/ 

r H 

2A 



(45) 



where Q < and / = ±2^-QA/3. 

Clearly, these cosmological solutions are explicitly dependent on r in a way which 
cannot be evaded by any coordinate transformation. As such they imply an intrin- 
sically inhomogeneous dark radiation efective energy density. Indeed, R is not a 
factorizable function, a fact preventing a reduction to the standard homogeneous dS 



or Robertson- Walker spaces @] 



3.2 AdS dynamics: A < 



In this scenario Eq. 



can be written as 



, 1 3 



dY 



2V A 7 

Now the only possibility is /3 > 0. Integrating we obtain the solutions 



^ +2A 



(3 sin 



/r2 -U 3/ ' 
-1 / ' 2A 



(46) 



(47) 



where for Q < the energy function can take any value in the range / > — 1 but for 



Q > it must also be in the interval |/| > 2J—QA/3. 



3.3 Absence of cosmological constant: A = 

In the limit A = Eq. (^) is written as 

RdR 



±t + ij = J 



Then the solutions for / 7^ are given by 



R' 



f 



Q+ ±ft + Jfr'-Q 



where fR^-Q > 0. Note that if Q > then / > 0. 



(4J 



(49) 



3.4 Absence of dark radiation: Q = 

When dark radiation is not present in the vaccum we have Q = and then the 
brane world dynamics depends on A and / as can be seen by going back to Eq. (^) . 
Integrating it we obtain 
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±t + i!{r) 



dR 



(50) 



Let us first consider A > 0. The solutions are now organized by /. For / > we find 



R 



If — 1 < / < we have 




sinh 



R = \ cosli 
V A 



±\ —t + sinh"^ 
V 3 



±J^t + cosh"^ 




(51) 



A 

3/' 



(52) 



Note that it is only for / = — Ar^/3 that solution (|52|) takes the homogeneous dS 
form [133 . 



When A < the relevant integral is only defined for / > and the corresponding 
solution is 



R 



-3/ . 
— : — sm 
A 



(53) 



In this case the homogeneous AdS space solution is again obtained with the choice 
/ = -ArV3. 



4 Physical Singularities and Regular Rebounces 

Both shell focusing singularities {R = 0) and shell crossing singularities {R' = 0) are 
found depending on the dynamics prescribed by Eq. (^91) . Shell crossing singularities 
are generally not considered to be "real" singularities, so we will concentrate always 
on shell focusing singularities in the following. 

To determine the inhomogeneous dynamical solutions corresponding to dark radi- 
ation gravitation on the brane we have implicitly assumed the validity of the following 
regularity condition 



V = V{R,r) = ^R'' + fR^-Q>0. 



(54) 



If Eq. (|54 ) is verified for all i? > then collapsing dark radiation shells will meet 
at the physical shell focusing singularity R = Rg = 0. This is what happens in the 
gravitational collapse of ordinary dust matter when A is set to zero . However, if 
A 7^ the evolution curve of a collapsing dust cloud has special rebouncing points at 
R > for which the dust potential is zero These points may exist both when a 
singularity forms at i?s = and when it does not. In the latter case the solution is 
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globally regular and the gravitational evolution of the dust cloud involves a collapsing 
phase followed by continuous expansion after velocity reversal at the rebounce point. 

Not surprisingly, an analogously rich structure of solutions can also be uncovered 
in the gravitational dynamics of vaccum dark radiation on the brane (see also the 
discussion in Ref. [0]). We start by writing Eq. (^) in the following way 

R^Bj" = V. (55) 

Starting from an initial (collapsing) state, a rebounce will occur for a given shell 
whenever = before that shell becomes singular. This happens when V{R,r) = 0. 
Introducing Y and /3 we write the potential V as the following simple quadratic 
polynomial 

V = ^(Y^-f3) (56) 



3 

and it becomes clear that no more than two regular rebounce epochs can be found. 
Depending on / and on the values of A and Q none, one or two such epochs may be 
associated with globally regular solutions. 

4.1 dS dynamics: A > 

The dS dynamics corresponding to A > is characterized by a concave phase space 
curve because 

d^V 2A , , 

This implies the existence of a phase of continuous expansion towards infinity with 
ever increasing speed. Other phases are possible depending on the parameters. 



For /3 < we have Q < and |/| < 2^-QA/3. Clearly, V is positive for all > 0. 
It satisfies V{0, r) = — Q and grows to infinity with R as Ai?^ (see Fig. m). The dark 
radiation shells may either expand continuously with ever increasing kinetic energy 
or collapse to a shell focusing singularity at Rs = after a proper time t = ts{r) given 

by 




sinh-i ( ^-^iM) - sinh-i ^ 



(5^ 



For A < Ac gravity is always confined near the brane during the cycle of evolution. 
However, if A > Ac then the gravitational field only stays confined for R < R^. As a 
consequence an expanding dark radiation vaccum goes through a phase transition at 
R = Rc where the gravitational field ceases to be confined to the brane. As we have 
seen this case and the large value A = Ac are currently ruled out by observations. 

For P > 0, Q > and / > — 1 (see Fig. ^ we find globally regular solutions with 
a single rebounce epoch at R = R^ where 
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Figure 1: Plots of F for A > 0, /? < and Q < 0. Non-zero values of / belong to the 
interval / > — 1 A I/I < 2^/— QA/S and correspond to shells of constant r. In the lefthand 
plot A < Ac and gravity is always bound to the brane. In the righthand plot A > Ac and 
the shaded region indicates where gravity is not confined near the brane. 



= + V^- (59) 

This is the minimum possible radius a collapsing dark radiation shell can have. The 
collapsing shells reach this epoch at the proper time t = t*(r) where 

At this point they reverse their motion and expand forever with ever increasing 
rate. The phase space of allowed dynamics defined by V and R is thus restricted to 
the region R > R^. Below R^ there is a forbidden region where the potential V is 
negative. In particular V{0, r) = —Q < which means that the singularity at i?s = 
does not form and so the solutions are globally regular. Note that if gravity is to stay 
confined to the vicinity of the brane for i? > i?* then A < Ac and R^, > R^. If not 
then there is a phase transition epoch R = R^ such that for R < R^ the gravitational 
field is no longer confined to the brane. For A > Ac gravity will always be free to 
propagate far away from the brane a situation that as we have seen is not allowed by 
current observations. 

If /3 > and Q < we have seen that / > — 1 has to further satisfy |/| > 



2y —QA/3. If / > 2y — QA/3 there no rebounce points in the allowed dynamical 
region R> and the evolution follows the case (3 < ilustrated in Fig. |1[ The time 
to reach the singularity at i?^ = is now given by 




, .... (^] 



cosh ^ — cosh ^ 



(61) 
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Figure 2: Plots of F for A > 0, /3 > and Q > 0. Non-zero values of / belong to the 
interval / > — 1 and correspond to shells of constant r. In the lefthand plot A > Ac gravity 
is always free to move away from the brane. In the righthand plot A < Ac and the shaded 
region indicates where gravity is not bound to the brane. 




On the other hand for — 1 < / < — 2y — QA/3 (see Fig. y) there are two rebounce 
epochs at i? = R^± with 

i?L = -^±V^- (62) 

Since V{0,r) = —Q > a singularity also forms at -R^ = 0. The region between 
the two rebounce points is forbidden because there V is negative. The phase space 
of allowed dynamics is thus divided in two disconnected regions separated by the 
forbidden interval i?*- < R < R*+- For < i? < the dark radiation shells may 
expand to a maximum radius R = R^^ in the time = where 



(63) 

At this rebounce epoch the shells start to fall to the singularity which is met after 
the proper time 

If i? > then there is a collapsing phase to the minimum radius R = taking the 
time = followed by reversal and subsequent accelerated continuous expansion. 
The singularity at -Rs = does not form and so the solutions are globally regular. 
Once more note that for A < Ac gravity is always confined near the brane but for 
A > Ac there is a phase transition epoch at R = Rc such that for R > R^. gravity is 
free to propagate out of the brane. Only if R^_ < R^ < R^+ this phase transition will 
not occur. For R < R,,_ gravity is confined near the brane and for R > it will 
not. 
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Figure 3: Plots of F for A > 0, /? > and Q < 0. Non-zero values of / belong to the 
interval — 1 < / < —2^J—QK/2> and correspond to shells of constant r. In the lefthand plot 
A < Ac and gravity is always confined near the brane. In the righthand plot A > Ac and 
the shaded region indicates where gravity is not bound to the brane. 




Figure 4: Plots of F for A > 0, /3 = and Q < 0. Here / = -2^-QK/'i. In the lefthand 
plot A < Ac and gravity is always confined near the brane. In the righthand plot A > Ac 
and the shaded region indicates where gravity is not bound to the brane. 



If ^ = then again Q < 0. For / = (see Fig. |) the single candidate 

to be a rebounce point is R = R^, with 

R* = (65) 

In this case ^(0, r) = —Q > and then a singularity also forms at Rs = 0. There is 
no forbidden region in phase space but the point at R^, turns out to be a regular fixed 
point which divides two distinct dynamical regions. Indeed if a shell starts at i? = i?* 
then it will not move for all times. If initially R < R^ then either the shell expands 
towards i?* or it collapses to the singularity. The time to meet the singularity is finite, 
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Figure 5: Plots of 1/ for A < 0, /? > and Q < 0. Non-zero values of / belong to the 
interval / > — 1 and correspond to shells of constant r. Here gravity is always bound to the 
vicinity of the brane. 




ts{r) 



2A 



In 



3Q 
A 



3Q 
A 



(66) 



but the time to expand to i?* is infinite. If initially R> then the collapsing dark 
radiation shells also take an infinite time to reach i?* 



If / = 2^-gA/3 there are no 
real rebounce epochs and the collapsing dark radiation simply falls to the singularity 
at Rs = 0. The colision proper time is 




(67) 



Of course for A < Ac then the gravitational field is always confined to the brane but 
for A > Ac there is a phase transition epoch R = R^ in either of the two dynamical 
regions such that for R > R^ gravity may freely move far into the bulk. 



4.2 AdS dynamics: A < 

The AdS dynamics corresponding to A < is characterized by a convex phase space 
curve because 

d^V 2A , , 

which does not allow a phase of continuous expansion to infinity. Also in this scenario 
the only possibility is /5 > 0. 

If Q < then we may consider / > — 1 (see Fig. Just like the corresponding 
dS case there is a single rebounce point at i? = i?* with 
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Figure 6: Plots of ^ for A < 0, /3 > and Q > 0. Non-zero values of / belong to 
the interval / > 2^/— QA/S and correspond to shells of constant r. In the lefthand plot 
Rc < R*- and gravity is confined near the brane in < R < R^^. In the righthand plot 
< Rc < R*+ and the shaded region indicates where gravity is not bound to the brane. 



Ri 



2A ^ V/^- 



(69) 



Due to the convexity of the phase space evolution curve this is the maximum radius 
which an expanding dark radiation shell can reach. The time for this to happen is 




TT 

sm 



V2 -U 3/ ' 
' 2A 



VP 



(70) 



Because V{0, r) = —Q > a singularity also forms at Rg = 0. Consequently, the 
reversal aX R = R^ will be followed by a collapsing phase to the singularity in the 
time 




7r 



— sm 



/ 3/ 



(71) 



2 V2Av^y. ■ 

This is in contrast with the dS scenario where there was no formation of a singularity 
at = and after reversal at the rebounce epoch the dark radiation shells expand 
forever. Also, because A < Ac gravity remains confined near the brane for all times. 



If Q > then / > 2y — QA/3 (see Fig. P). We find two regular rebouncing epochs 
at -R = i?*± where 



Rl± 



2A 



(72) 



The allowed phase space region is precisely the interval between these two roots, 
R^,- < R < R*+. Since the singularity at i?s = is inside the forbidden region it does 
not form and so we are in the presence of globally regular solutions. In this case an 
initially expanding dark radiation shell will rebounce at R = -R*+ after the time 
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71 

2 - 



(73) 



Then it reverses its motion to collapse until it reaches R = in the time 




(74) 



It rebounces again to expand to i? = in the same amount of time and subse- 
quently repeats the cycle. Thus these are oscillating globally regular solutions. Again 
note the contrast with the dS scenario where no oscillating solutions were found. 
Since A < Ac there is only a gravity confinement phase transition epoch in the cycle 
if < Rc < 



4.3 Absence of cosmological constant: A = 

When A = the potential is simply V = fR^ - Q. If Q < and / > then there 
are no rebounce epochs in the allowed dynamical range R> and the singularity is 
at i?5 = where V{0,r) — —Q > 0. The collapsing dark radiation shells reach the 
singularity after the time 

ts{r) = J {Jfr^ - g - v^) . (75) 
In this case gravity never propagates out of the brane. If Q > and / > or Q < 



and — 1 < / < there is a single rebounce root at i?* = \JQ/ f- For Q > and 
/ > the allowed dynamical region is R> R^ and so the singularity does not form 
at i?5 = 0. A collapsing shell falls to a minimum radius in the time 

U{r) = J (76) 

and then expands continuously to infinity. If R^ > Rc gravity is confined for all 
R > Rj,. Otherwise the dark radiation goes through a confinement phase transition 
at R = Rc- For Q < and — 1 < / < the dynamical range in phase space is 
< R < R^. The shells may expand out to reach R^ at 



= J ■ (77) 

Then it reverses its motion to fall into the singular epoch i?^ = in the proper time 

tsir) = (78) 
In this setting gravity remains confined in the allowed dynamical range, < R < R^. 
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4.4 Absence of dark radiation: Q = 

When the dark radiation degrees of freedom are not present Q — and the potential 
to be considered is 



V = V{R,r) = ^R^ + f 



(79) 



For A > we may distinguish three types of evolution. If — 1 < / < there is 
a rebounce epoch towards infinite continuous expansion at -R* = 



'-3// A and the 

allowed dynamical region is > i?*. The time to collapse to minimum radius is 



tJr) 



-cosh ^ 



■37. 



(80) 



If / > there are no rebounce epochs for > and a singularity forms at Rg — 
where the potential is V{0, r) — f. The time to meet the singularity is 



tsir) 



-sinh ^ 




(81) 



If A < the only possibility is to take / > 0. Then we find a rebounce epoch at 



R* 



3//A. After expansion to a maximum radius i?* in the time 



3 
A 



TT 



— sm 



-1 



r\ -- 



A ' 



(82) 



the brane collapses to the singularity at Rs — reaching it after the time 



TT 6 

tJr) = -\ --. 



A 

For A < Ac the gravitational field is always confined. Otherwise it never is. 



(83) 



5 Conclusions 

We have analysed some aspects of the gravitational dynamics of inhomogeneous dark 
radiation on a RS brane. This is an important first step towards the understanding 
of a realistic collapse setting on a brane. Indeed, the behaviour of dark radiation is 
nothing else than the dynamics of the vaccum on the brane when the gravitational 
modes are excited only as energy density. 

We have taken an effective 4-dimensional view point to show that with certain 
simplifying, but natural assumptions, the Einstein field equations form a solvable, 
closed system. The solutions obtained were shown to depend on the brane cosmo- 
logical constant A, the dark radiation tidal charge, Q, and on the energy function 
fir). 
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We have given a precise description of the dynamics of the solutions and charac- 
terized how this depends on A, Q and f{r). We have also presented the conditions 
defining the solutions as singular or as globally regular. 

Finally we have discussed the confinement of gravity to the vicinity of the brane. 
We have seen that it depends on A, Q and also on the brane tension A. If A < Ac 
where Ac = and Q < then the gravitational field is always confined to the 

vicinity of the brane during its evolution. Alternatively, if A < Ac and Q > gravity 
is only bound to the brane for R > where i?^ = 3Q/(Ac — A). For A > Ac and 
Q < confinement is restricted to the epochs R < Rc- We have noted that current 
observations do not allow A < and A > Ac. Consequently, the gravitational degrees 
of freedom may eventually propagate away from the brane at high energies. Such a 
phase transition would occur at the epoch R = R^ and should be detected in a change 
in Newton's gravitational constant. 
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